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i ' Abstract 

<: 

r \ \ There is a contrast between the two sets of functional equations 

>S : fo(x + y) = fo(x)My) + fi(x)h(y), 

fi(x + y) = f 1 (x)My) + f (x)f 1 (y), 

l — '■ and 

^ ! fo(x -y) = fo(x)f (y) - /i(x)/i(y), 

>; h(x-y) = h(x)fo(y)- fo(x)h(y) 

satished by the even and odd components of a solution of f(x + y) = f(x)f(y). 
J. Schwaiger and, later, W. Forg-Rob and J. Schwaiger considered the extension 
of these ideas to the case where / is sum of n components. Here we shorten and 

q^ \ simplify the statements and proofs of some of these results by a more systematic 

0^\ ' use of matrix notation. 
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1 Introduction 

By a more systematic use of matrix notation, we shorten the statements and proofs of 
some results of J. Schwaiger ^H] o n generalized hyperbolic functions and their char- 
acterization by functional equations. We also discuss some stability results for these 
equations; our results here are motivated by those of W. Forg-Rob and J. Schwaiger 
PI; see also [TUj . 

The results discussed here can be readily understood in the case n — 2. A function 
/ : C —* C can be written as a sum of even and odd components in a standard way: 
fix) = f (x) + A{x), where f (x) = [f(x) + f(-x)}/2, f ± (x) = [f(x) - f(-x)}/2. In 
case / is an exponential function, i.e, if / satisfies 

f(x + y)=f(x)f(y), (1.1) 

for all x and y, then the components satisfy 

fo(x + y) = fo(x)f (y) + fi(x)fi(y), n 9 ^ 

h(x + y) = A(x)/ (y) + f (x)h(y), { ^ } 

and 

fo(x -y) = fo(x)fo(y) - fi(x)fi(y), . . 

h{x -y) = h(x)fo(y) - fo(x)h(y). { 6> 

In case / is the usual exponential function, then f and f\ are the cosh and sinh 
functions, respectively, and ([1.2)1 . ([1.3)1 are the familiar sum and difference relations 
for these functions. But, of course, ([1.1)1 also has many discontinuous solutions [TJ 
Chapters 2 and 3]. 

Conversely, the general solution of (jl. 3[) is expressible in terms of a single arbitrary 
exponential function, i.e., it is known [20] that if f and /i satisfy ([1.3)1 . then f 
and /i are the even and odd components of a single exponential function: fo(x) = 
[g(x) + g(-x)}/2, fx(x) = [g(x) - g(-x)]/2. In fact, g(x) = f (x) + fx(x). On the 



other hand |19| . the general solution of (jl-2|) depends on two exponential functions: 
fo(x) = [gi{x) +g 2 (x)]/2, /i(x) = [g-i (x) - g 2 {x)]/2. We could think of g x (x) as e x , for 
example, and g 2 as the identically function. 

The reason that (|1.3jl has fewer solutions than (J1.2J1 is that ()1.3|) implies (J1-2J1 but not 
vice- versa. To see that (|1.3|) implies (jl.2|) . we first interchange a; and £/ in (|1.3jl to see 
that /o must be even and /i odd. Then replacing y by — y in (jl.3j) we get (jl.2j) . On 
the other hand, /o(;e) = /i(x) = 6^/2 satisfies (J1.2J1 but not (|1.3|) . 

We get "higher order generalized hyperbolic functions" by considering a function writ- 
ten as a sum of n components, forming a natural extension of even and odd components 
in the case n = 2. The present work discusses the appropriate generalizations of (jl.2|) . 
(J1.3J) and their solutions in this setting. Theorem 2.8 gives the systems of equations 
analogous to (J1.2J1 . (jl.3j) and Theorem 3.1 gives the general solutions of these systems. 
Theorem 5.2 discusses the stability of these systems and makes it clear why the system 
(jl.3|) is stable but the system (jl.2j) is not. 

The name "higher order generalized hyperbolic functions" is intended to convey that 
the functions discussed here generalise the hyperbolic functions in two ways. First, 
there is the general order n rather than the special order 2. Second, we do not confine 
attention to continuous solutions of the relevant functional equations. 



The main results reported in §2 and §3 are proved in [To], essentially by the same 
methods, but with a less systematic use of matrix methods. An attempt to understand 
[TH] and [Oj is what led to the present paper. In (THj and [Oj, n-tuples of functions 
were used; here we find it more convenient to base the development more explicitly on 
circulant matrices formed from these n-tuples; this is partially done in ^Bj and [H]- 
The question of stability is further considered in ^Hl ■ 



We note that (jl.2j) and ()1.3|) can be written 



and 



fo(x + y) fi{x + y) 
fi(x + y) fo(x + y) 



fo(x-y) fi{x-y) 
fi(x-y) fo(x-y) 



My) h{y) 
h(y) fo(y) 



fo{y) -h{y) 
-fi(y) fo(y) 



fo{x) fi(x) 
fi(x) fo(x) 



fo(x) f x {x) 
fi(x) fo(x) 



;i-4) 



;i.5) 



so our task is to find higher-dimensional analogues of (J1.4JI and f)1.5|) . 



Our results are stated for complex- valued functions of a complex variable, but, as in [9, 
[TU] and PHI, they can be extended to complex- valued functions on an abelian group 
G with an automorphism a satisfying a n = ida- 

F. Zorzitto [21] has considered a far-reaching generalization, based on group theory, of 
the results of |H] and [TH] . 

2 Preliminary results 
2.1 Special matrices 



The n x n Fourier matrix [7J p. 32] T n is given by 



J~ J~n. 



UJ 



[0 



CO 



1 UJ-( n -V OJ-^-V 



u -in-l) 
w -2(n-l) 



U 



-{n-lf 



Tn is a symmetric matrix whose complex-conjugate is given by 



J~ — J~ n 



1 1 

1 

1 



10 

uj 2 



iV 4 



1 L) n - 1 LU 2( - n -V 



u n - 1 

w 2(n~l) 



U> 



(n-1) 2 



(2.1) 



Sometimes the Fourier matrix is denned as in (j2.1[) but without the factor 1/y/n. The 
present definition has the advantage that the matrix T n is unitary: T^Tn = Fn^n = I- 
We use the notation |7J circ(a T ) = circ(ao, a±, . . . , o n _i) for the circulant matrix: 



A = circ(a T ) = circ(a , ai, . . . , a n _ 



ciq a\ ... a n _i 

On-l O-O ■ ■ ■ 0"n-2 
CL\ 0,2 ■ ■ ■ ^0 



(2.2) 



We also use the diagonal matrix 



Q = Q n = diagfl, u, . . . , uj 



n-ln 



We have Theorem 3.2.1] 



F*VLT = 7T, 



(2.3) 



or 



Q = TtxT* 



(2.4) 



where n is the n x n permutation matrix 



7r = 7r n = circ(0, 1, . . . , 0) 



In fact all circulants are polynomials in 7r; we have [TJ p. 68] 
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irc(o ,ai,. . .,o n _i) =p 1 (n), 



cue 



(2.5) 



where 



p-y{z) = a + aiz H h a n _i2; 



n-l 



(2.6) 



Every circulant matrix can be diagonalized: 



Lemma 2.1. If A is a circulant matrix given by \2.ty) . then 

TAT* = v^ diag[.F* a]. (2.7) 

Here we use the usual boldface notation for column vectors: a = [a a x ... a n _i] T 
and diagfa] denotes the square matrix whose diagonal elements are do, a%, . . . , a n -\- 

Proof. This is shown in ([7, Theorem 3.2.2]) in the form 

TAT* = diag(p 7 (l),p 7 (u;), . . . ,p 7 (^" 1 )), (2.8) 

which is equivalent to (J2.7J) on using (j2.fi j) . 

By conjugation and iteration of ()2.3|) we get, since n T = n^ 1 and T T = T, 

TQ- m T* = n m , (2.9) 

for m = 0, 1, . . . and, taking inverses, this also holds for m = —1, —2, .... This gives 

T[n~ rn 7i k Q m ]T* = TQ- m T*[Tn k T*}Tn m T* = n m Tn k T*n~ m , (2.10) 

for integers m and nonnegative integers k. Hence, in view of the representation ()2.5|) . 
we get: 

Lemma 2.2. For circulant A, and m — 0, 1, . . ., 

T[n- m AQ m }T* = 7i m TAT*n~ m . (2.11) 



2.2 Components of a function 

Throughout the rest of this paper, we let n be a fixed integer > 2. Following 
J. Schwaiger [TH] . we make the following definition: 



Definition 2.3. A function h : C — > C is of type j (j = 0, . . . ,n — 1), if h(ujx) = 
ujihiyx) where uj = uj n = e 2m l n . 

In case n = 2, the type and type 1 functions are the even and odd functions, respec- 
tively. In the cases n = 3 and n = 4, examples of functions of types 0, 1, 2 and 0, 1, 2, 3, 
respectively, are given in (|4.1jl and (J4.2J) . It is easy to generate examples in the case 
of analytic functions: if n — 3 and f(z) = Y^kLo a kZ k i in a disk centred at 0, then 
fj( z ) — J2kLo a 3k+j z3h+:j is of type j, j = 0, 1, 2. But we have no need of analyticity in 
the discussion. 



Lemma 2.4. Every function f : C — >• C can be expressed uniquely as a sum of 
functions fj of type j, j — 0, . . . , n — 1, called the components of f . In fact, 

71-1 



where 



fo(x) 
fi(x) 

fn-l(x) 



T 



n 



3=0 



f(ux) 



f{uo n - l x) 



T* 



n 



f\u^x) 
f{oj~ n+1 x) 



(2.12) 



;2.i3) 



Proof. Equation (J2.12J) with fj of type j, j — 0, . . . , n — 1, implies that 



/H 

f^u^x) 



^r 



fo(x) 
fi(x) 

fn-l(x) 



[2.14) 



which is obtained from (|2.12jl by replacing x by x, ux, . . . , c<j n_1 x in order. Also, of 
course, (J2.14J) implies (J2.12J) . But (|2.14|) obviously has a unique solution for the fa 
given by the first equality in (j2.13|) . and each of the fi is of type i. 



The second equation in (|2.13|) is obtained similarly. 



Definition 2.5. The circulant matrix 



/o /l • • • /n-l 

/n-1 /O • • • /n-2 



(2.15) 



F(s) =circ(/o,/i,...,/„-i) 

/i h ■ ■ ■ /o 

whose first row is formed from the components of f in order will be called the circulant 
matrix corresponding to /. 

Lemma 2.6. /// : C — ► C zs any function, the corresponding circulant matrix function 
is given by 

F(x) = F*diag [/(*), /H, . . . , /(w- 1 !)]^ (2.16) 

Proof. This comes from Lemma 2.1. 

2.3 Exponential functions 

Lemma 2.7. /// is an exponential function, the corresponding circulant matrix func- 
tion F satisfies F(x + y) — F(y)F(x). 

Proof. This is a fairly immediate consequence of Lemma 2.6. J-F(x + w)^ 7 * is a 
diagonal matrix with entries f(uj l (x + y)) = f( y u t x)f(uj' t y), i = 0, ...,n — 1. But 
f(cu l x),f(u! l y),i = 0,...,n — 1, are the respective entries of the diagonal matrices 
FF(x)F* and FF(y)F*. 

Lemma 2.8. For an?/ function f , the corresponding circulant matrix function F sat- 
isfies the equation F(uy) = Q~ 1 F(y)Q. 

Proof. This is a fairly immediate consequence of Definition 2.3. 

Successive application of Lemmas 2.7 and 2.8 give 

Theorem 2.9. If f is an exponential function, the corresponding circulant matrix 



function F satisfies 

F{x + to m y) = n~ m F(y)n m F(x), (2.17) 



for m = 0, 1, . . . , n — 1. 



Theorem 2.9 is a matrix form of jTHJ equation^)], which is written in the notation 

3 n-l 

f j (x + cu t y) = '£ujW t f l (x)f M (y)+ £ u^'* fi(x)f „*-&), t = 0, 1, . . . ,n - 1. 
z=o J=j+i 

In case n = 2, f)2.17|) reduces to ()1.4|) for m = and to (|1.5jl for m — 1. 

3 A converse result 

We now consider equation (j2.17|) where -F(x) is circulant, but no assumption is made 
about its entries /j. It will come as no surprise that if (|2.17|) holds for all m = 
0,1,. . . ,n — 1, then the /j are the components of an exponential function and F(x) = 
circ(/o, /i, . . • , fn-i)- But what if we have (J2.17J) for only a single value of m? Then it is 
clear that (J2.17J) holds also for integer multiples of this value. Thus if g.c.d.(n, m) = 1, 
and (|2.17|1 holds for a fixed m 6 {0, 1, . . . , n— 1}, then it holds for all m G {0, 1, . . . , n — 
1}. On the other hand, if (|2.17|) holds for a fixed m G {0, 1, ... ,n — 1} with g.c.d. 
(n,m) = d > 1, then we can assert only that it holds for values of m which are 
multiples of d, modulo n, and it turns out that the general solution of ()2.17|) depends 
on d arbitrary exponential functions: 

Theorem 3.1. Let a circulant matrix function F(x) satisfy 

F(x + u m y) = tt- m F(y)n m F(x), (3.1) 



for a fixed value of m G {0, 1, . . . , n — 1}, and let d = g.c.d.(n, m). Then F(x) 
circ(fYx) T ) where 



fix) = -^T 



n 



h(x) 

h(u; m a;) 



(3.2) 



h ^m(n/d-l) x ^ 

h(x) being an arbitrary d-tuple of exponential functions. In particular, when d = 1, 



fix) = -=T 



9{x) 

g(ujx) 



g(uj n l x) 



where g is an arbitrary exponential function. 



To prove Theorem 3.1, we use the ideas of [EJ. We write (J3.1J) in the form 



(3.3) 



FF{x + u m y)F* = [m- m F{y)n m T*][TF{x)T* 



(3.4) 



and using Lemmas 2.1 and 2.2, this is 



G(x + uo m y) = n m G(y)7r- m G(x) 



(3.5) 



where G(x) is the diagonal matrix 



G(j/) = ^[diag^*f(j/)]. 



Using the notation 



G m (x) = Tt m G{x)Tj- 



(3.6) 



we have 



G m (x) = dia,g[g m (x), g m+ x(x), . . . , g m + n -i{x)], 



where suffixes are taken modulo n, since 



(3.7) 



7rdiag[Ai, . . . X n }iT = diag[A 2 , . . . A„, A 



i > 



so the effect of the transformation G(x) — > 7i m G(x)n m is to replace Qj by gj. 
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The easiest case is m = and hence d = g.c.d.(m, n) = n. In this case (13. 5|) is 

G(x + y) = G(x)G(y), (3.8) 

so, since G is diagonal, its entries must be exponential functions. In case m = this is 
all we can say about the g's; they are independent. 

The other extreme case is where d — 1. In this case the set of multiples {km, k = 
0, 1, . . . , n — 1} modulo n is precisely the set {1, 2, . . . , n} so ()3.5|) can be written as a 

system of n equations 

g (x + u m y) = g m (y)g (x), 
g m (x + uj m y) = g2m(y)g m (x), / 3 9 n 

3(r.-i)m(a: + u m y) = g (y)9(n-i)m(x). 
There are two possibilities. One is that go(0) = 0; hence from ([3.9)1 . all of the <& are 

identically and the assertion of the Theorem holds trivially. The other possibility is 

that <7o(0) 7^ 0. Putting x = y = in the first of equations (|3.9|) then gives g m (0) = 1. 

Then the remaining equations, with y = 0, give <?2m(0) = <?3m(0) = • • • = go(0) = 1 or 

#2(0) = #3(0) = • • • = go(0) = 1. Then putting x = in (J3.9J) we see that the g% are 

given by the formula 

tffcmG/) = 9o(oo km y), k = 0, 1, . . . , n - 1. (3.10) 

The second of these equations gives 

g m (y) = 9o^ m y), (3.n) 

which, together with the first of equations ()3.9|) shows that go is an exponential function. 
Now (|3.1U|) is equivalent to 

g k {y)=go(u k y), k = 0,1,..., n- 1. (3.12) 

Thus G is determined by the single exponential function g and f is determined by 
()3.3|) . This completes the proof in the case m — 1. 
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In case d = g.c.d.(n, m) > 1, the equations (J3.9J) are replaced by d sets of n/d equations 
corresponding to the cosets of the additive subgroup of {0, 1, ... ,n — 1} generated 
by {0, 1, . . . , d — 1}. The rth such set of equations (where r is one of the numbers 

0, 1, . . . , d— 1), reads 

g r (x + u m y) = g r+m (y)g r (x), 
g r+m (x + u m y) = g r+2m (y)g r + m (x), , g ^ 

g r +(n/d-l)m(x + UJ m y) = g r +(n/d)m(y)gr+(n/d-l)m(x) ■ 

Each such set of equations can be dealt with as in the case m — 1. That is either all 
of the g r +km (k — 0,1, ... , n/d — 1) are identically or they are given by the formula 



g r +km(x) = g r (uj m x), k = 0, 1, . . . , n/d 



(3.14) 



where g r is an exponential function. Thus G is determined by the ci-tuple h(x) T = 
(g (x) , gi(x) , . . . ,g d _i(x)) of exponential functions and we get the main assertion of 
Theorem 3.1. 

4 Continuous solutions 



In the case n — 1, a continuous solution of f(x + y) = f(x)f(y), is given by f(x) = e ax . 
For n > 1, e x has components 

-r n 

^2(*) = -E^' fc ^, fc = 0,... l n-l, 

in a minor variation of the notation of [T7j . |18j . For further references and history 
of these functions, see JH] • It should be noted that the term "generalized hyperbolic 
functions" as used in [TJ], [T2j and [IB], corresponds to some of the continuous solutions 
of the system of functional equations considered in the present paper. The three 



continuous generalized hyperbolic functions of order 3 are 



*So(s) = 


1 
3 


& x + 2e -x/2 cog v|r" 


? 


3,1 \ X ) = 3 


e x_ 2e -s/2 C0S (v|z + !y 


3,2( 2 ') = 3 


e x 


- 2e- x / 2 cos (^ - 


1)1 



(4.1) 
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In the case n = 4, we get the compact formulas 

F} Q {x) = (1/2) (cosh x + cos x), 
Fli( x ) = (l/2)(sinhx + sinx), 
Fl 2 (x) = (1/2) (cosh a; — cosx), 
Fl 3 (x) = (l/2)(sinhx — sinx), 

given by Battioni [4^ , a special case of a more general result 

FLA X ) = l F kr( x ) + F ™U x )}/ 2 , r = 0,l,...,m-l, 

^2m,r+m\ X ) = [^ m,r\ X ) ~ ^ m,r\ X )\l ^^ 7" = U, 1, . . . , m — 1, 

given in a different notation in (33), p. 216]. 

The decomposition of various special functions into sums of functions of order j, j = 
0, . . . , n — 1 has been considered by Y. Ben Cheikh; see j^j and [§]. 



5 Stability 

We use the usual 1-norm for square matrices: 

n 

PII = Plli= ma? XX- 1, 
If B is a diagonal matrix, we have 

||AB|| > ||A||min 1 <j< n |6 M |. (5.1) 

Definition 5.1. (See, e.g., |3], |I3j-) For fixed m, we say that the equation 

F(x + u m y) = Q- m F(y)n m F(x), (5.2) 

for circulant F , is stable if the condition 

\\F{x + u m y) - n- m F(y)tt m F(x)\\ < e, for some e > and all x, y, (5.3) 

implies that either F is bounded or F satisfies \5.ty) . 
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It has been noted that in the case n = 2, the system (|1.3|) is stable JT3J but (jl.2|) is not 
. An appropriate generalization is 



Theorem 5.2. Let F be a circulant n x n matrix function and let m be a fixed integer 
in the set {0, 1, . . . ,n — 1}. The equation 

F(x + u m y) = Q~ rn F{y)Q m F{x), (5.4) 

is stable if and only if g.c.d.(n,m) = 1. 

Proof. The idea of the proof is based on the proof of Proposition 2]. Suppose that 
g.c.d.(n,m) = 1 and let F be unbounded and satisfy (|5.3|) . We need to show that F 
satisfies ()5.4|) . We introduce the diagonal matrix 

G(x) = TF[x)T* = Vn diag [JF*f ( x )\ 

as in Lemma 2.1. Clearly G is also unbounded, since ||JF|| = ||.F*|| = ^/n. These 
relations, together with Lemma 2.2, show that (|5.3jl is equivalent to 

\\G(x + u m y) - G m (y)G(x) \\<ne = e 1 , (5.5) 

with the notation of (|3.6j) . (J3.7J) . This leads to 

\\G(x + (u m y + u m z)) - G m (y + z)G{x) || < e x 

and 

||G((x + u m z) + u m y) - G m (y)G(x + u m z) \\ < e x 

so 

\\G m (y + z)G(x) - G m (y)G(x + co m z) \\ < 2e 1 , 

for all x,y,z. Also, from (|5.5p and the norm property ||^4-B|| < ||A||||S||, we get 

\\G m (y)G(x + u m z)-G m (y)G m (z)G{x)\\<e 1 \\G r 
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Combining the last two inequalities, we get 

\\G m {y + z)G(x)-G m (y)G m (z)G{x)\\<e 1 (2 + \\G m (y)\\). 

Since G is diagonal, we may use (|5.1jl to get 

\\G m {y + z) - G m (y)G m (z)\\mm.i<i< n \gi(x)\ < \\[G m (y + z) - G m (y)G m (z)]G(x)\\ 

< s 1 (2+\\G m {y)\\). (5.6) 

We see from ()5.5jl that \\G{uj m y) — G m (y)G(0)\\ < e\\ hence the unboundedness of one 
diagonal entry in G implies the unboundedness of all diagonal entries in G. (Here it 
is crucial that g.c.d.(n, m) = 1; otherwise integers multiples of m will not cover all 
j = 0, . . . ,n — 1.) This means that, given y and z, we can choose x such that 

min 1 < i < n |p i (a;)| > 2+ ||G m (?/)||. 

Hence, from ()5.6|) . \\G m {y + z) — G m (y)G m (z)\\ < e\. Since the elements of G m are the 
same as those of G, in a different order, we finally get 

\\G(y + z) - G(y)G(z)\\ < e 1 . (5.7) 

Thus each element gi in the diagonal matrix G(y) is unbounded and satisfies the in- 
equality \gi(y + z) — gi(y)gi(z)\ < £\. By [21 Theorem 1], we must therefore have 
9i(y + z)- gi(y)gi(z) = 0, so 

G(y + z)-G(y)G(z) = (5.8) 

and, as in the proof of Theorem 2.9, F is a solution of (J5.2J) . 

To prove the "only if assertion of the Theorem, we note that if g.c.d.(n, m) > 1, 
the unboundedness of one diagonal entry in G, does not imply the unboundedness 
of all entries there. So, for example, suppose that gi(x) = e x , for i a multiple of 
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d = g.c.d. (n,m) > 1 and gi(x) = 2, when % is not a multiple of d. Then G is an 
unbounded solution of (|5.7jl . with £1 = 2, which does not satisfy ()5.8|) . Thus we can 
have an unbounded solution of (J5.3J) which does not satisfy (J5.4J) . This, of course, is 
exemplified by the non-stability of the equations (jl.2|) . 

Acknowledgments. I thank the referees for very helpful suggestions and for identi- 
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